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Abstract 

We investigate both the hyperbolic action and the determinant ratio action designed to fix 
the topological charge on the lattice. We show to what extent topology is fixed depending on 
the parameters of these actions, keeping the physical situation fixed. At the same time the 
agreement between different definitions of topological charge - the field theoretic and the index 
definition - is directly correlated to the degree topology is fixed. Moreover, it turns out that the 
two definitions agree very well. We also study finite volume effects arising in the static potential 
and related quantities due to topology fixing. 



1 Introduction and motivation 



Phenomena related to topology of gauge fields are of inherent non-perturbative nature. Using 
lattice techniques to investigate topology is, therefore, a natural choice. However, a lattice 
discretization leads to fundamental problems, e.g. naively discretizing the topological charge 
density yields non-integer values for the topological charge in general. Moreover, with different 
discretizations significantly different results can be encountered rendering it difficult to establish 
a unique concept of topology on the lattice. 

To stress the importance of topology, we mention the mass of the rl meson, which originates to 
a large extent from topological effects. In addition, the low lying Dirac eigenvalues are sensitive 
to the topological charge [1], which is also visible in random matrix theory [2] and in chiral 
perturbation theory in the so-called e-regime [31 SI [5]. 

Therefore, it is highly desirable to have a well defined notion of topology on the lattice. One 
approach in this direction is to use the standard continuum expression of the topological charge 
density, q oc FF ^ with a lattice discretization of the field strength, in particular improved versions 
of it. This is combined with suitable smearing techniques smoothing the sampled gauge field 
configurations a 'posteriori. Another approach is to use the index of the overlap operator, which 
leads to a conceptually clean integer definition by means of the index theorem (but may also 
depend on parameters of the kernel of the overlap operator). While it has been demonstrated 
that smearing and other filtering methods agree to a large extent on the locations of topological 
structures [6l|7l[8], the agreement of the total topological charge is not satisfactory for precision 
measurements. 

A different approach to arrive at a unique definition of topological charge is to restrict the gen- 
eration of lattice gauge configurations a priori to those, which are smooth enough to resemble 
continuous gauge configurations. It has first been shown by Liischer p] that the topologi- 
cal charge has an integer valued so-called geometric definition on lattice configurations, whose 
quantum fiuctuations are bounded, i.e. whose plaquettes are close to the identity. In order to 
implement such a constraint on the plaquette values, modified gauge actions can be used such 
as the hyperbolic action of ref. |10] . It has been demonstrated in practical simulations that the 
use of the hyperbolic action (and variants of it) indeed tends to fix topological charge |1H I12|. 

Constraints on the plaquette that guarantee the locality of the overlap operator |131 ITi] represent 
similar admissibility conditions. Obeying these conditions assures that the kernel of the overlap 
operator has a spectral gap and hence the eigenvalue density of zero or near-zero eigenvalues 
vanishes. As a consequence, zero mode crossings and in turn changes of topological charge are 
avoided. 

Such a method has been implemented in jl5| with the aim to perform simulations with dynam- 
ical overlap fermions. The idea is to introduce a determinant ratio employing Dirac operators 
at negative quark mass [16| . which leads to the desired spectral gap and hence forbids topology 
changes. This is, because current dynamical overlap algorithms including approximations of 
the quark determinant make changes in the topological sector (i.e. overcoming the correspond- 
ing discontinuity) very costly, which results in a scaling with 0{V'^) |17j . When performing 
simulations with fixed topological charge, there is no need to treat topological tunneling and 
simulations with dynamical overlap fermions become feasible as they are currently carried out 
by the JLQCD collaboration. For a recent review of these simulations see [18], where it is also 



1 



demonstrated that topology is indeed fixed and that the method outhned above can be used in 
practice. 

While the simulations of the JLQCD collaboration aim at obtaining QCD results, which could 
be compared to experiment, we want to address a more theoretical question. It is our goal to 
investigate, whether using a topology fixing action (either the hyperbolic action or the deter- 
minant ratio action) will lead to an agreement regarding the values of the topological charge, 
when different lattice definitions are used. Besides the theoretical and conceptual interest in this 
question a practically very useful consequence could be that much less computer time demand- 
ing definitions of the topological charge than the index of the overlap operator can be used. 
This would open the possibility to address topology on much larger lattices than it is presently 
feasible with the overlap operator. A first investigation concerning the agreement of different 
definitions of topological charge using the hyperbolic action has been performed in 

In this work, we will concentrate on two definitions of topological charge, the field theoretic 
definition via an improved field strength tensor, which also involves some additional smearing, 
and the fermionic definition via the index of the overlap operator. As we will show below we 
find that using topology fixing actions goes hand in hand with a better agreement between these 
two lattice definitions of topological charge. 

It is important to note that working at a fixed value of the total topological charge in a finite 
volume violates cluster decomposition. In sufficiently large volumes local topological fluctuations 
still take place in every fixed topological sector. Therefore, fixing topology on the lattice to e.g. 
the trivial sector is supposed to provide correct physical results in the infinite volume limit with 
cluster decomposition becoming restored. We expect that fixing topology induces particular 
finite volume effects in accordance with [191 I20j . where these effects have been addressed by 
using the ^-dependence of physical quantities. 

Due to these finite volume effects physical observables receive different values, when evaluated in 
different topological sectors in finite (and small) volumes. In this paper we address the question 
by a preliminary investigation of the dependence of Wilson loops, the static potential and the 
Sommer parameter tq. We indeed find clear evidence that such a dependence is present. We 
stress that it is the original idea and long term plan of the present project to adress topological 
finite size effects in a detailed manner computing many more observables than the ones studied 
here. We consider this paper as a first and essential step in this direction and we believe that it is 
an important basic work for such a future investigation. Our main and novel result that for the 
determinant ratio action different topological charge definitions agree, allows to determine the 
topological charge also on very large lattices where the overlap definition of topological charge 
becomes impractical. 

Our results indicate that the value of the topological charge will not depend on the details of 
its definition in contrast to findings in the past on rather coarse configurations, without any 
smoothing. We find this result quite remarkable and important. The field theoretic definition is 
known to be strongly affected by local lattice artefacts in general (and how they are treated e.g. 
under smearing) . With the determinant ratio action we employ a method that constrains a global 
observable (the lowest eigenvalue of the kernel of the overlap operator) almost by construction. 
However, as we show for the first time here, it also leads to unique values of topological charges 
through the field theoretic definition. 

The paper is organized as follows. In the next section we write down two definitions of topological 
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charge on the lattice. The two sections afterwards are devoted to the implementation of two 
topology fixing actions, the hyperbolic action and the determinant ratio action. We compare 
obtained results to those from the standard Wilson action. In section [S] we investigate the 
dependence of the static potential and related quantities on the topological sector. We close 
with a brief comparison and conclusions. 



2 Definitions of the topological charge 
2.1 Field theoretic definition 

The straightforward definition of topological charge on the lattice is the discretization of the 
standard continuum expression, 

= ^EWI^(^M<^(^)^P<x(x)), (1) 

X 

where F^y{x) is some lattice field strength tensor consisting of appropriate linear combinations 
of Wilson loops. Instead of the naive discretization of the field strength tensor, which combines 
only plaquettes, i.e. the smallest possible Wilson loops we employ an improved field 

strength tensor [21] 

F,,{x) = hCl^J<'\x) + k2C;ix'\x) + hCjix'\x), (2) 

where C^"^"'*(x) are the clover averages of n x n Wilson loops attached to the lattice site x. 
In order to achieve O(a^) improvement at the tree level, we choose ki = 1.5, k2 = —0.15 and 
k3 = 1/90. 

In addition, smoothing or filtering techniques are needed to come closer to integer values for Q, 
i.e. reducing renormalization factors originating from short-range fiuctuations. We use 10 steps 
of APE smearing before evaluating eq. ([TJ. APE smearing \23\ [M] is an iterative procedure, 
where links are replaced by weighted averages of the old links C/^ and their attached staples 
= Uy{x)U^{x + i')ul{x + fi) projected back to the gauge group: 

U, ^ Psu(3){{l-a)U^ + ^^U;'\. (3) 

The projection onto SU{3) is not unique. Psu{3)(^) projects W onto that V € SU{3) that max- 
imizes Re{Tv{VW^)), where the maximum is found iteratively as implemented in the CHROMA 
software for lattice QCD [25]. Setting the weight parameter a = 0.45, APE smearing has been 
argued to be equivalent to RG cycling [26] . 

Smearing suppresses UV fluctuations in every observable, but is also biased towards classical 
solutions in the long run. Even in the short run it may destroy small topological objects. Hence, 
10 steps can be considered as a compromise between obtaining topological charges closer to 
integer values and the destruction of topological objects. 

^ Stout smearing has been shown to lead to very similar results for the topological charge [H]. 
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2.2 Definition via the index 



An alternative and conceptually clean way of defining the topological charge motivated by the 
index theorem in the continuum is to use the index of the overlap operator [27^ I28j . The 
overlap operator, which is a solution of the Ginsparg- Wilson relation [29], has been introduced 
by Neuberger [301 [3T] and reads for massless quarks 

L>ov = 1 + , A = Dw-mo, (4) 

where I?w is the Wilson Dirac operator and we have chosen the mass parameter ttiq = 1.6 (as 
in |15| 132)). In the massless limit, which we consider here, the overlap operator (cf. eq. Q) has 
exact zero modes with a definite chirality. This allows to define the topological charge as 

Q' = ^^Tr(75Z)ov(x,x)) = index(Z)ov) = N_ - N+ , (5) 

X 

where -^-(+) denotes the number of zero modes with negative (positive) chirality. To identify 
the zero modes we follow the techniques described in 



2.3 Comparison of topological charge definitions on standard Wilson gauge 
configurations 

Let us give an example for the difference of the resulting topological charges, when using the two 
definitions given above. We consider gauge configurations generated with the standard Wilson 
action 

Sw[U] = f3Y,Sp{Up) , Sp{Up) = l-Up , Up = iRe(Tr(t^(i,i))), (6) 
where Up denotes the plaquette variable and /3 = 6/gQ. 

In the left panel of fig. [1] we plot the Monte Carlo time histor}{^ of the difference of the two 
topological charge definitions, 

AQfc = Qi-Ql (7) 

as function of the configuration index k using /? = 6.18 (this corresponds to a ~ 0.07 fm) and 
a lattice volume of 16^. We also show the corresponding histogram in the right panel of that 
figure. As can clearly be seen frequent non-zero values of AQk are encountered. This shows that 
the two definitions of the topological charge yield different results. In particular, configurations 
with |AQfc| > 0.5 would by rounding be ascribed to different topological sectors. Interestingly, 
we observed that most of these configurations are associated with transitions in the index along 
the Monte Carlo time and that during such transitions the field theoretical definition typically 
recognizes the new topological charge not as quickly as the index definition does. 



One unit of Monte Carlo time corresponds to ten sweeps of link updating. 
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Figure 1: The difference of the topological charge AQk = Qk~Qk as function of the configuration 
index k (left) and the corresponding histogram (right) (Wilson plaquette gauge action). 

Since we are already considering highly improved definitions of the topological charge, we do 
not follow the line of improving the topological charge or the filtering even further. 

An alternative route might be the application of improved actions. For example with the concept 
of perfect actions [35] it can be expected that different definitions of the topological charge lead 
to coinciding values. 

Another possibility, which we will follow in this work, is to use actions that fix topology. Such 
actions have been tailored to address questions related to topology of gauge fields and give rise 
to the expectation that again better agreement of different topological charge definitions will be 
reached. Particular realizations of topology fixing actions are the hyperbolic action discussed 
in [361 [T2] and actions employing a determinant ratio of Wilson- type fermions [ll^ [TS] . We will 
give the precise form of these actions in the next two sections. We will also discuss, whether 
using these actions indeed leads to a better agreement of the topological charge obtained from 
different definitions. 



3 The hyperbolic action 

The hyperbolic action provides a cut-off for small plaquette values (i.e. large action density). It 
is given by 



^-^ ' ' ' ' ' oo otherwise 



where for the rest of the paper we will use n = 1. This form of the action has been suggested and 
used in [371 EB] for conceptual studies of chiral gauge theories on the lattice and it is motivated 
by the locality studies of the overlap operator [13]. There it has also been shown that for 
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the Wilson kernel of the overlap operator has a spectral gap and hence zero mode crossings 
and, therefore, changes of topology are forbidden. Note that the bound in eq. Q has later been 
improved [U] to 

Sp(Up) < e = ^ ~ (10) 

^ ^ (1 + l/y2)d(d - 1) 20.5 ^ ^ 

As we will describe below these theoretical bounds are typically not reached in practical simu- 
lations and, therefore, the fixing of topology will be imperfect. 

Note that the form of the action (cf. eq. ([8])) does not allow the application of a heatbath 
algorithm. Therefore, it has been simulated by using a local Hybrid Monte Carlo algorithm as 
proposed in [39] and discussed in its application to the present case in |12| . 

The aim of our simulations is to use the hyperbolic action (and the determinant ratio action 
to be discussed later) in the same physical situation (i.e. same physical volume and lattice 
spacing) defined by the example of the Wilson gauge action discussed in section [2] and fig. [TJ 
Consequently, when decreasing the value of e, also the value of /3 has to be decreased in order 
to stay at the same physical lattice spacing as was demonstrated in |12j . 

In fact we took the same combinations of parameters (/3, 1/e) detailed in [12], but performed 
independent simulations. The results for the plaquette expectation value (Up), the integrated 
plaquette autocorrelation time T^^^'^ and the acceptance rate and the step size of our local HMC 
algorithm are given in Table [TJ In addition we list the frequency of topology changes /top in the 
corresponding runs, where we measured the topological charge through the index of the overlap 
operator. In the last column we give the static force parameter ro/a [4^ as computed from 
smeared Wilson loops using standard procedures [41 j . From this the lattice spacing in physical 
units can be inferred, using ro = 0.5 fm. 





lA 


(Up) 


plaq 
'int 


acceptance rate 


dT 


/top 


ro/a 


6.18 


0.0 


0.61185(2) 


8.2(9) 


0.998 


0.1 


0.0223 


7.17(9) 


1.5 


1.0 


0.599858(5) 


2.2(2) 


0.997 


0.1 


0.0066 


7.09(8) 


1.0 


1.18 


0.601518(3) 


1.20(6) 


0.995 


0.1 


0.0047 


7.30(7) 


0.8 


1.25 


0.598366(3) 


1.13(5) 


0.992 


0.1 


0.0038 


7.08(8) 


0.3 


1.52 


0.601035(2) 


0.82(2) 


0.877 


0.1 


0.0008 


7.5(2) 



Table 1: Simulation parameters and results for the hyperbolic action (eq. dH]), lattice volume 
16^). (Up) denotes the plaquette expectation value, t^^^^ the corresponding integrated plaquette 
autocorrelation time, dr the Monte Carlo step size, /top the frequency of topology changes and 
ro/a the static force parameter. 

Table [T] clearly shows that, when staying in the fixed physical situation of ro/a ~ 7, which 
amounts to a lattice spacing of approximately 0.07 fm, the frequency of topology flips is getting 
smaller, when decreasing e as has already been noted in [12]. Note, however, that even at our 
smallest value of e = 0.66 (1/e = 1.52 in Tabled]) topology is not completely fixed. 

In fact at ro/a ~ 7 it will not be possible to fix topology completely. The qualitative argument 
for this is the following: a roughly constant plaquette expectation value can be viewed as an 
indicator of being in a fixed physical situation. This is also suggested by Table [1) a value of 
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ro/a ^ 7 corresponds to {Up) ^ 0.6. Decreasing e beyond the desired value (Sp) is impossible, 
since according to eq. ([8]) the weight of all configurations with action density Sp > e becomes 
zero. Hence, there is a critical value ecrit ~ (•S'p), which is the lowest value that can be reached for 
a given physical situation. For our present case e cannot be lowered below ecrit ~ 1 — 0.6 = 0.4. 

In Table [2] it is demonstrated that going to larger values of rg/a (i.e. to smaller values of the 
lattice spacing) the frequency of topology changes is reduced as expected. This indicates that at 
sufficiently small lattice spacings the hyperbolic action will lead to a complete fixing of topology. 
However, this will presumably happen at so small values of the lattice spacing that numerical 
simulations will be extremely demanding. 
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lA 


dT 


/top 


ro/a 


0.25 


1.52 


0.01 


0.0041 


6.7(2) 


0.3 


1.52 


0.01 


0.0012 


7.5(1) 


0.1 


1.64 


0.01 


0.0048 


6.4(1) 


0.15 


1.64 


0.01 


0.0007 


8.1(2) 



Table 2: Simulation parameters and results for the hyperbolic action of eq. ([S]). Going to smaller 
values of the lattice spacing (i.e. to larger values of rg/a) amounts to decreasing the frequency 
of topology changes /top. 

As a consequence of the above discussion our situation of ro/a ~ 7 can only be achieved, when 
using values of e>0.4. In practice, we stopped our simulation at e = 1/1.52 0.66, because 
for smaller values it is rather difficult to do a fine tuning of e and/or /? to the desired physical 
situation. Naturally the question arises, whether for this value of e there is already agreement 
between the topological charge values obtained from the two different definitions explained in 
section [2j As measure for the (dis)agreement we use 

where the difference AQ^ is defined in eq. d?]). 

In the left panel of fig. [2] we show the frequency of topology changes and in the right panel of 
the same figure the values of AQ as functions of /?. AQ exhibits a similar behavior as ftop- it 
decreases as topology is fixed more and more. We conclude that topology fixing actions have 
indeed the potential to assure a unique results for the topological charge from different defini- 
tions. Remember, however, that at typical values of the lattice spacing it is not possible to fix 
topology sufficiently strongly with the hyperbolic action. Therefore, we consider an alternative 
approach using the determinant ratio action. 

4 The determinant ratio action 

The idea of fixing topology originates from the locality property of the overlap operator |13j . 
For a sufficiently fine lattice spacing (a < 0.1 fm) the kernel operator has a spectral gap and 
hence the overlap operator is local. This can be proven analytically using the theoretical bounds 
on e for very small values of the lattice spacing. It has also been checked numerically. 



7 



0.03 
0.025 

0.02 
0.015 

0.01 
0.005 




rn/a = 7 



3 4 



O 



0.1 
0.08 
0.06 
0.04 
0.02 




•i 



1 2 3 4 5 6 7 



Figure 2: The frequency of topology changes /top (left) and the difference of the topological 
charge AQ (right) as functions of j3 for approximately fixed r^/a (hyperbolic action). Note 
that we do not show statistical errors for /top, because of the long autocorrelation time of that 
quantity; nevertheless we consider the value of /top as an indicator to what extent topology is 
fixed. 



A spectral gap of the kernel operator in turn implies that there are no zero mode crossings and 
hence a change of topology (when using the index definition) is forbidden. This line of arguments 
suggests to include in the partition function a term 

det (a^A) 

Rdct = (12) 

det (^A^A + j 

Note that here the kernel operator A (cf. eq. ([4])) is evaluated at negative quark mass —mo to 
guarantee the existence of a spectral gap. In the denominator in eq. (jl2p a regulating twisted 
mass term /.ttm has been introduced, which controls the range of near zero modes that are to be 
suppressed. 

The determinant ratio action has been suggested in |1HI15]. There it has also been checked that 
it indeed fixes topology, when the topological charge is determined by the index of the overlap 
operator. The determinant ratio action has been extensively used in simulations with dynamical 
overlap fermions (see |18] for a recent review) to avoid slowing down the algorithm at changes 
of the topological charge (see section [T]). 

Here, we use the determinant ratio action as a local modification of the Wilson plaquette action. 
We do not include any dynamical fermions and, therefore, work in the quenched approximation. 
Our main motivation is to compare the field theoretic and the index definition of the topological 
charge and to determine to what extent they yield consistent results. 

The determinants that appear in eq. (jl2p are treated by means of pseudo fermionic fields and a 
global Hybrid Monte Carlo algorithm |42j . Details of this HMC algorithm can be found in |43] . 
It uses even/odd preconditioning 01], mass shift (l5lll6] and multiple time step integrators. 

In a first step we have confirmed the results in [111 [T5] that the topological charge, when 
computed via the index of the overlap operator, is indeed fixed. To this end we generated initial 
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configurations on 16^ lattices with topological charge = 0,1,2 and performed simulations 
with the determinant ratio action. We have chosen /3 = 6.063, mo = — 1.6 and ^tm = 0.2 (see 
below). We found that on samples of 300 configurations separated by 10 trajectories not a single 
topology change took place. Thus, topology seems to be fixed. 

To be able to perform a meaningful comparison with the hyperbolic action results from section [3l 
we tuned the value of /3 such that ro/a ~ 7. The results of this tuning are collected in Table [3l 
where we list /3, the average plaquette (Up), the plaquette integrated autocorrelation time Tj^'f^, 
the acceptance rate and the force parameter ro/o. Note that these results correspond to the 
= sector. The main results of this table are that a value of /3 = 6.063 leads to the desired 
scale of ro/a ~ 7 and that the change in f3 to reach this physical situation is much milder than 
in the case of the hyperbolic action. 





(Up) 


plaq 
'int 


acc. rate 


ro/a 


5.900 


0.58850(3) 


3.7(6) 


0.85 


5.43(4) 


6.000 


0.59877(3) 


3.0(5) 


0.89 


6.40(9) 


6.050 


0.60358(3) 


2.9(5) 


0.90 


6.69(6) 


6.063 


0.60482(2) 


2.3(3) 


0.90 


6.97(5) 


6.100 


0.60830(3) 


2.5(4) 


0.91 


8.04(22) 



Table 3: Simulation parameters and results for the determinant ratio action of eq. ()12p . 

The crucial question is, of course, whether the two definitions of the topological charge, (cf. 
eq. (dJ) and (cf. eq. agree. In fig. [3] the /3 dependence of AQ (eq. (fTT]) ) is shown for 
the simulations detailed in Table El Note that in all these runs = 0, i.e. all simulations 
were performed in the topologically trivial sector. One can clearly see that AQ decreases, when 
approaching the continuum limit by increasing f3. For (3 = 6.063, i.e. for ro/a ~ 7, AQ ~ 0.02 
indicates rather good agreement between the two definitions. 

0.1 
0.08 
0.06 

o 

<i 

0.04 
0.02 


5.85 5.9 5.95 6 6.05 6.1 6.15 
P 

Figure 3: AQ as function of /? for the determinant ratio action. Note that contrary to fig. [2] 
here we change the physical situation, i.e. the value of tq/u (cf. Table [3]). 

In the upper panels of fig. [J] we show the Monte Carlo history of AQi^, the disagreement of the 
two definitions, at /? = 6.063. When rounding the field theoretic results to the nearest integer, 
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there is agreement for almost all gauge configurations. Only for the lower values of (3 we have 
encountered in rare cases differences bigger than 0.5. Corresponding plots for the hyperbolic 
action (/3 = 0.3, 1/e = 1.52) look essentially identical (cf. the lower panels of fig. S]). 




Figure 4: The difference of the topological charge AQ^ = Q\^—Q\ as function of the configuration 
index k (left) and the corresponding histogram (right) for the determinant ratio action at /? = 
6.063 (upper panels) and for the hyperbolic action at /3 = 0.3, 1/e = 1.52 (lower panels). 

We close this section with a short remark about autocorrelation times. As can be seen in 
Table [3] the autocorrelation times of the plaquette turn out to be rather short in our simulations. 
However, it is important to realize that the plaquette does not provide the largest autocorrelation 
time of the system. This can be seen from fig. [5l where we show the autocorrelation times of 
smeared Wilson loops (extension R x T) as functions of R keeping T = 7 fixed. Note that the 
autocorrelation time in fig. [5l tint , is measured in units of configurations and each configuration 
is separated by ten Monte Carlo iterations. Therefore the autocorrelation times of the Wilson 
loops in fig. [5] are an order of magnitude larger than the ones for the plaquette given in Table [3l 
Clearly, autocorrelation times become larger, when increasing R until at ~ 6 they level off 
and reach plateaus. The increase of the smeared Wilson loops autocorrelation times is not a 
property characteristic for the topology fixing actions only [47] . In particular, when computing 
the static potential and from that ro/a one should use these larger autocorrelation times, when 
computing statistical errors. 
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Figure 5: Autocorrelation times of smeared Wilson loops with extension R functions of 

R with T = 7 fixed, here shown for a fixed topological sector, Q = (determinant ratio action, 
/3 = 6.063, lattice volume 16^). Note that the autocorrelation time in fig. [5] is measured in units 
of configurations and each configuration is separated by ten Monte Carlo iterations. 



5 Dependence of the static potential and related quantities on 
the topological sector 

Another important issue, which we intend to address in more detail in a subsequent publication, 
is the dependence of physical observables on topology, in particular the finite size effects in a 
given topological sector. Such results could be confronted with the theoretical expectations of 

Here we study in a first step the Q dependence of the static potential and related quantities. We 
employed the determinant ratio action at /3 = 6.063 using 16^ lattices. To identify the topological 
charge we used the index definition Q^. We computed smeared Wilson loop averages (VI^(_r,t)) 
(extension R x T), the static potential V{R) at separation R and the Sommer parameter ro 
in the topological sectors with charges S {0 , 1 , 2 , 4}. As can be seen from Table IH there 
are substantial differences between results computed in different topological sectors going well 
beyond a 5a effect, in particular when comparing = and = 4. This is also visible from 
fig. El while the static potentials for different sectors look rather similar (left), a zoomed plot 
shows a clear dependence on the topological charge (right). 

In a recent paper by the JLQCD Collaboration |32| a similar analysis was carried out: tq 
was computed for the topological sectors S {0,2,4}, but no dependence on within 
statistical errors has been established. Note, however, that compared to our simulations the 
spatial extension of the lattice L (in physical units) was larger by a factor of about 1.7 and 
the temporal extension larger by a factor of ~ 3.4, yielding in total a factor of about 17 for 
the lattice volume V = L"^ x T. Assuming that the dependence of tq is of the same form 
as that derived for spectral quantities in \19\ |20] the absence of a Q-dependence is consistent, 
because the dominating finite volume correction, when fixing topology, is proportional to 1/V, 
which presumably would shrink the deviations seen here to the 1 to 2a level. Note also that 
in |32] dynamical quarks were used, which, in principle, could have an effect since topological 
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(W^(4,4)> 


(W^(7,7)) 


y(4) 


V{7) 


ro/a 





0.94776(2) 


0.4090(6) 


0.0931(7) 


0.2177(3) 


0.3326(8) 


6.97(5) 


1 


0.94770(3) 


0.4074(6) 


0.0911(8) 


0.2189(3) 


0.3364(9) 


6.84(4) 


2 


0.94765(2) 


0.4057(6) 


0.0895(7) 


0.2201(3) 


0.3393(9) 


6.77(4) 


4 


0.94743(2) 


0.4000(4) 


0.0832(5) 


0.2245(2) 


0.3515(7) 


6.48(3) 



Table 4: Dependence of various observables related to the static potential on the topological 
sector defined by Q^. (VF(_r,t)) denotes smeared Wilson loop averages (extension R x T), V{R) 
the static potential at separation R and tq the Sommer parameter (determinant ratio action, 
/3 = 6.063, lattice volume 16^). 
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Figure 6: The static potential as a function of the separation for the topological sectors = 
{0 , 1,2,4} (left), and a zoomed plot of the same results (right). 

fluctuations are suppressed. Moreover, the lattice spacing used in this reference was significantly 
larger (0.12 fm compared to 0.07 fm). 



6 Discussion and summary 

In this paper we were mainly concerned with the question, whether different definitions of 
the topological charge lead to consistent results, when using topology fixing gauge actions. 
To address this question, we employed two modifications of the Wilson plaquette action, the 
hyperbolic action, which implements Liischer's admissibility condition, and the determinant 
ratio action, which leads to a spectral gap. 

As a result we could establish that there is a correlation between the degree topology is fixed 
and the agreement of both topological charge definitions used herein in the sense that the better 
topology is fixed the better the agreement comes out. Let us accept for a second the index of the 
overlap operator at a given value of tuq as a reference definition of the topological charge, which 
is not unreasonable given the conceptual cleanness. Then we find that the hyperbolic action is 
not able to fix topology completely, even for the smallest value of e ~ 0.66 reached in this paper. 
We argued that in fact with a value of ro/a ~ 7 it will not be possible to achieve such a complete 
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fixing of the topological charge. On the other hand the determinant ratio action fixes topology 
at once for ro/a>5.4, i.e. already at rather coarse lattice spacings. Note that in [32] it has been 
found that even at ro/a ~ 4.2 no change of the topological charge could be observed. Of course, 
it is still possible that modifications of the overlap kernel operator such as changing the value of 
niQ or smearing the gauge configurations may lead to different values of the topological charge 
(the index), a question that clearly deserves further studies. 

When fixing topology in a finite volume, as is necessary for certain numerical simulations, a 
number of conceptual questions arise. As a first step towards addressing such problems, we 
studied the dependence of the static potential on the topological charge. For our lattice size 
with linear extent of L ~ 1.2 fm we found a clear dependence of the static potential at distances 
of about R ~ 0.5 fm, when computed in the topological charge sectors Q = and Q = 4. It will 
be very interesting to see, whether this difference vanishes when taking the infinite volume limit. 
It is our goal to investigate the volume dependence of physical observables in fixed topological 
sectors. This will also allow us to test the predictions for the finite size effects given in \19\ I20j. 

In this work we focused on global aspects of topology. It would also be interesting to investigate, 
how topology is locally realized in topology fixed gauge configurations. Moreover, we expect 
that the mechanism of topology fixing is rather different for the two actions we have used. This 
can be seen from fig. \7\ where we show the plaquette distributions for the Wilson, the hyperbolic 
and the determinant ratio action at the same physical situation. 

As shown in more detail in fig. 3 of |12j the hyperbolic action shows an asymmetric distribution 
for l/e> 1.52, since small plaquette values are suppressed. In fact this suppression has been the 
original motivation for fixing topology. In contrast to that the distribution of the determinant 
ratio action is very similar to the distribution of the Wilson plaquette action. This can be 
made more quantitative by fitting Gaussian functions to these distributions with mean values 
/i and widths a as parameters. The results of these fits are listed in Table [5] confirming that 
the value of a is almost smaller by a factor of two for the hyperbolic action. Nevertheless, the 
determinant ratio action fixes the reference topological charge completely. This result suggests 
that the suppression of small plaquette values is not a necessary condition for fixing topology. 
It would be interesting to investigate this in more detail. 



action 




lA 




a 


^hyp 


0.8 


1.25 


0.598360(9) 


0.000254(4) 


Sdet 


6.063 




0.604808(9) 


0.000479(7) 


Sw 


6.18 




0.611857(9) 


0.000497(7) 



Table 5: Results from Gaussian fits oc exp a*)^/2o-^ plaquette distributions shown in 

fig.El 

Acknowledgments 

We thank Gunnar Bali, Gregorio Herdoiza, Michael Miiller-Preussker, Silvia Necco, Luigi Scorzato 
and Andreas Schafer for useful discussions and in particular Urs Wenger for helpful cross checks 
of static potential computations. 

This work has been supported in part by the DFG Sonderforschungsbereich/Transregio SFB/TR9- 



13 



600 


P=0.8 


|3=6.063 


500 


p=D. 1 O 
1 


400 




300 




200 




100 










0.595 0.6 
Hyperbolic 



0.605 0.61 0.615 

Detratio Wilson 



Figure 7: The plaquette distributions for the Wilson, the hyperbohc (/3 = 0.8) and the deter- 
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